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We theoretically investigate the effect of broken inversion symmetry on the generation and con¬ 
trol of ultrafast currents in a transparent dielectric (Si 02 ) by strong femto-second optical laser 
pulses. Ab-initio simulations based on time-dependent density functional theory predict ultrafast 
DC currents that can be viewed as a non-linear photogalvanic effect. Most surprisingly, the direc¬ 
tion of the current undergoes a sudden reversal above a critical threshold value of laser intensity 
Ic ~ 3.8 X 10^^ W/cm^. We trace this switching to the transition from non-linear polarization cur¬ 
rents to the tunneling excitation regime. We demonstrate control of the ultrafast currents by the 
time delay between two laser pulses. We find the ultrafast current control by the non-linear pho¬ 
togalvanic effect to be remarkably robust and insensitive to laser-pulse shape and carrier-envelope 
phase. 


In the last decade, ultrafast few-cycle laser pulses with 
well-defined carrier-envelope phase (CEP) have become 
available providing novel opportunities to explore the ul¬ 
trafast and non-linear response matter to strong optical 
fields. The study of the induced electronic motion and 
of the highly non-linear optical response have focussed 
on rare gas atoms [l[, molecules [H and, more recently, 
on nanostructures, surfaces and bulk matter Hi- The 
driven electron dynamics can be monitored through op¬ 
tical signals and through emitted electrons [lol-llii . 
Very recently, Schiffrin et al [l6[ have demonstrated di¬ 
rected electron currents generated inside transparent di¬ 
electrics by carefully tailored laser pulses. In turn, the ul¬ 
trafast response can characterize the impinging laser field 
la. Currently, avenues are explored to exploit such ul¬ 
trafast modulation of electric currents for petahertz-scale 
signal processing [isl enabled by the short intrinsic time 
scale of the electron motion (^ 1 fs), orders of magnitude 
faster than semiconductor electronics. 

In this letter, we explore a novel channel for the ul¬ 
trafast electronic response that is unique to dielectrics 
with a non-centrosymmetric crystallographic structure: 
the generation of DC currents induced by strong opti¬ 
cal laser pulses. Fully three-dimensional ab-initio simu¬ 
lations based on time-dependent density functional the¬ 
ory (TD-DFT) predict the generation of strongly non¬ 
linear currents in a-quartz that are, in contrast to pre¬ 
viously observed currents independent of the 

details of the laser pulse shape. The direction of the 
currents is found sensitive to the instantaneous laser in¬ 
tensity. Analysis of the spatio-temporal charge dynamics 
on the atomic length and time scale allow us to link this 
to the transition from non-linear polarization currents to 
directional tunneling excitation, the latter being highly 
sensitive to the alignment between the laser polarization 
and the chemical bonds in the crystal. We demonstrate 


that this transition may be investigated in a pump-probe 
setup leaving its marks as a change of the direction of the 
current as a function of the pump-probe delay. 

Theoretical exploration of ultrafast processes in solids 
faces the challenge to tackle the time-dependent many- 
body problem. Time-dependent density functional the¬ 
ory has emerged as a versatile tool allowing for an ab- 
initio description of a variety of strong field processes 
in the solid state [1, 13, 15, 0. Here, we ernploy a 
real-space, real-time formulation of TDDFT |2QI425| for 
the electronic dynamics induced by strong few-cycle laser 
pulses in (a-Si02 (a-quartz). Briefly, we solve the time- 
dependent Kohn-Sham equations (atomic units are used 
unless stated otherwise) 


idti^i{Y,t) = H{r,t)'ipi{r,t) , 


( 1 ) 


where i runs over the occupied Kohn-Sham orbitals ipi. 
The Hamiltonian 

Hir,t) = ^{-iV + A{t)f+Vion+j dr'^l^+rxc(r,i) 

( 2 ) 

describes the system under the influence of a homogenous 
time-dependent electric field F(t) of amplitude Fq along 
a with vector potential A(t) = — F(f)df in the ve¬ 
locity gauge and in the transverse geometry allowing 
to treat the bulk polarization response of the infinitely 
extended system along the polarization direction. The 
periodic lattice potential Von is given by norm-conserving 
pseudopotentials of the Troullier-Martins form repre¬ 
senting the ionic cores (O(ls^) and Si(ls^2s^2p^)). The 
valence electron density is n(r,t) = 
the exchange and correlation potential Vxc we employ 
the adiabatic Tran-Blaha modified Becke-Johnson meta- 
GCA functional [28|. It accurately reproduces the band 
gap A ^ 9 eV for Si02 and yields good agreement 







2 


with the experimental dielectric function over the spec¬ 
tral range of interest including at optical frequencies [29| . 
The time-dependent Kohn-Sham equations (Eq. [1]) are 
solved on a Cartesian grid with discretization ^ 0.25 
a.u. in laser polarization direction and ^ 0.45 a.u. per¬ 
pendicular to the polarization direction in a cuboid cell 
of dimensions 9.28x16.05x10.21 a.u.^ employing a nine- 
point stencil for the kinetic energy operator and a Bloch- 
momentum grid of 4^ k-points. The time evolution is 
performed with a 4^^-order Taylor approximation to the 
Hamiltonian with a time step of 0.02 a.u. including a 
predictor-corrector step. The solution of Eqs. [T] and [2] 
allows to analyze the time and space dependent micro¬ 
scopic vectorial current density 

= \e\Y^^[ip*{r,t){-iV + A{t))ipi{r,t) + c.c.] 

i 

(3) 

as well as the mean current density J{t) along the laser 
polarization direction Fq, averaged over the unit cell of 
volume 

J{t) = ^ J^dr -Fo/lFol . (4) 

The polarization density P{t) = gives 

the charge density D(t) transferred by the pulse. The 
total charge Q will depend also on the details of the ge¬ 
ometry of the laser focus and of the collection volume not 
explicitly treated in the following. 

Eirst studies of the short-pulse induced current and 
charge transfer in polycrystalline Si 02 EEllii found 
a sinusoidal dependence on the carrier-envelope phase, 
0CE5 of the few-cycle electromagnetic field A(t) ^ 
Ao(t) cos(ci;Lt-h 0 ce) with (cos^) envelope Ao(t) and ujl 
the carrier frequency of the IR laser. Subcycle control 
and steering of electrons required exquisite control over 
the instantaneous electric field F(t). Here, we explore an 
alternate route towards steering, controlling, and switch¬ 
ing ultrafast currents that does not rely on (/)ce control of 
the instantaneous field but on the instantaneous intensity 
dependence of a direct (DC) current. 

Starting point is the observation that the total charge 
density D{rp) transferred at the conclusion of the pulse 
can be split into a CEP dependent part with amplitude 
Dcep{^p) and a residual part Do{rp). Dcep{^p) tends 
to decrease with increasing pulse length while the mag¬ 
nitude of Do{rp) increases with the pulse length (Eig.[T]). 
Eor pulse length exceeding a few optical cycles, Do{rp) 
is approximately proportional to the pulse duration and 
dominates the signal exceeding Dqep by about one order 
of magnitude. 

The charge transferred by the induced DC current, 
Do(rp), features a strongly non-linear scaling with in¬ 
tensity I Do I (X or, equivalently, field strength |Do| oc 
(Eig. [2j). The origin of this highly non-linear re¬ 
sponse lies in the broken centrosymmetry of the Si 02 
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FIG. 1. (Color online) Pulse length dependence of transferred 
charge at intensities (a) 1 x 10^^ W/cm^, (b) 5 x 10^^ W/cm^, 
and (c) 1 X 10^^ W/cm^, each with cos^ pulse shape. Carrier- 
envelope phase (CEP) dependent part Dcep (blue open 
squares); CEP independent part Dq (green full squares); 
dashed lines: linear slope going through the origin. 



FIG. 2. (Color online) Carrier-envelope phase independent 
transferred charge density Do (green full squares) as function 
of laser intensity (cos^ pulse with full duration Tp — 20 fs, 
hwL = 1.7 eV) and absolute value |Do| (empty boxes). Power 
law |Do| oc (dashed line). Electrons move along -ha for 
I < Ic — 3.8 X lO^^W/cm^ (upper panel) and along —d for 
I ^ Ic (broken line to guide the eye). 


crystal along the d direction. In general, generation of a 
directed flow of charge by a laser field requires a broken 
inversion symmetry. Eor few-cycle laser pulses with well- 
defined CEP, inversion symmetry is violated by a suit¬ 
able choice of 0 ce- In the present case, it is no longer 
the temporal shape of the laser electric field but the elec¬ 
tronic and crystallographic structure of matter the laser 
interacts with that causes ultrafast currents. This novel 
mechanism does not rely on delicate CEP control yet of¬ 
fers sub-cycle response and switching. 

The appearance of a direct current in a homogeneous 
medium under illumination, independent of the CEP, and 
linearly increasing with pulse duration, can be viewed 
as a non-linear analogue to the well-known photogal- 
vanic (PG_) effect |3ll435j as first discussed qualitatively 
by Alon [36|. Conventionally, the lowest order photogal- 
vanic effect is described as 

= HklnFlF* . (5) 

JPG 

is quadratic in the electric field components and 
linear in the time-averaged laser intensity I oc FiFp 
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FIG. 3. (Color online) (a) Time-dependent polarization den¬ 
sity P{t) (solid green), laser field F{t) (black dashed, scaled) 
for Tp = 20 fs, ujl = 1.7 eV, intensity 1 x lO^^W/cm^. 
(b) Time-dependent non-linear current AJnl (Eq. [6l red 
solid), (c) a-quartz lattice with broken inversion symmetry 
(a —a). “Larger” atoms are closer to the reader than 
small atoms; labeled atoms form a helix along a direction, 
(d-e) Snapshot of the current density in an a-c cut plane 
going through the central oxygen at a time Q (d), @ (e), 
dash-dotted rectangle in (c). Angles between electric field 
and 0-Si bond are 71 = 25.3 deg, 72 = 51.5 deg. (f) Time- 
averaged density modulation after the laser pulse An(Tp) in 
plane cut through the oxygen marked by a star in (c) for 
intensity 1 x lO^^W/cm^. 


For linearly polarized light, the photogalvanic tensor 
l3kin associated with the two-wave mixing in the second- 
order susceptibility Xkini^l^ nonzero only in non- 

centrosymmetric crystals [37| . Microscopically, a variety 
of mechanisms may contribute to the PG effect such as 
asymmetric excitation, scattering, or recombination of 
electrons and electronic defects [32|. One important re¬ 
alization is the so-called “shift current” 
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due to 

the shift between the center of charge of the valence elec¬ 
trons and the excited electrons in the conduction band. 
This shift current has been predicted to be important in 


several semi-conductors [37. ,|4Q| and has been first exper¬ 
imentally verified for ferroelectrics [s^. 

The present non-linear generalization of photogalvanic 


effects is obviously a signature of strong-held interaction 
with matter. This is underscored by the surprising ob¬ 
servation of current reversal as a function of laser inten¬ 
sity (Fig. [2]). We hnd a critical value of current reversal 
at Ic = 3.8 X lO^^W/cm^. Electrons move along the 
+d direction for low intensities I < Ic while they prop¬ 
agate along —a direction for higher intensities I > Ic. 
We elucidate the microscopic mechanism for this rever¬ 
sal by analysis of the spatio-temporal charge dynamics. 
At lower intensities / < Ic, the multi-photon driven non¬ 
linear polarization currents lead to a localized accumu¬ 
lation of charge in between the Si-0 bonds as displayed 
in the time-averaged density huctuations at the conclu¬ 
sion of the laser pulse (Fig. [3f). This implies the forma¬ 
tion of an induced atomic-scale dipole around the oxygen 
atoms, i.e. a displacement of the center of charge by ver¬ 
tical excitation, resembling the shift current mechanism 
of the standard photogalvanic effect but generalized to 
higher order reflected in the non-linear intensity scaling 
of \Do\ ^ For higher laser intensities I > P the 

dominant charge transfer mechanism is excitation of the 
tilted conduction band by tunneling. Tunneling signifi¬ 
cantly depends on the local potential landscape in tun¬ 
neling direction. We find tunneling is enhanced when the 
bond direction is aligned with the laser field as evidenced 
by a strongly asymmetric current density at times near 
the maxima of the electric field (Fig.[3li,e). Tunneling ex¬ 
citation is more efficient along — d where the 0-Si bond is 
more closely aligned with the laser field (71 = 25.3 deg) 
while in +d direction tunneling is suppressed because of 
the larger angle (72 = 51.5 deg). The transition to tun¬ 
neling excitation is therefore accompanied by a reversal 
of the charge transfer and current direction. As the tun¬ 
neling rate scales exponentially with the peak intensity 
(X exp(—1/\/7), the transition is quite abrupt suggesting 
its potential for femtosecond current switching. 

For high intensities / > Ic, the charge transfer shows 
sub-cycle time structure. The time-dependence of the 
tunneling current can be conveniently analyzed by the 
non-linear response contribution AJnl(^) after subtract¬ 
ing the linear-response current scaled to the instanta¬ 
neous field, 

AJnl(^) = ./(i) - — / e “V(a;)F(w) (6) 

with the conductivity a{uj) determined for low intensity 
I < Ic [HI. During the rise time of the pulse (Fig. 
AJnl is still ^ 0 as linear response prevails. However, 
once a field strength is reached sufficient for tunneling be¬ 
tween neighboring atoms at around —3 fs, the non-linear 
current shows strong spikes. While the linear response 
current is, to a good approximation, 90 degrees out of 
phase with the electric field and P{t) = is 

in phase with F(t), the current spikes are in phase with 
F{t) as expected for tunneling excitation. At later times 
(from -1 fs), AJnl remains in phase with, but becomes 
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FIG. 4. (Color online) (a) Laser-induced transferred charge 
density as a function of the peak-peak delay {h = 2.4 x 
10^^ W/cm^, h = 0.6 X 10^^ W/cm^, ujl = 1.7eV, Tp = 20 fs). 
Red squares, dashed line: TD-DFT simulations; solid line: 
model Eq. [3 Inset: magnified data around At ~ 0 to 2fs 
(dash-dotted arrow), (b) Maximum instantaneous laser in¬ 
tensity I{t) — F(t)^(c/87r) for a given delay. Insets: pulse 
shape after superposition of pump and probe pulse F{t) at 
At = — 13fs (weak pulse before strong pulse) and at +10 fs. 


proportional to the laser field, indicative of a conductor¬ 
like linear response J(t) = (7DF{t) with a Drude (free 
carrier-like) conductivity an for the tunneling-induced 
electron population in the conduction band. 

The present analysis of the non-linear photogalvanic 
DC current suggests that the key control parameter is 
the instantaneous intensity /(t) rather than the cycle av¬ 
eraged intensity in the conventional photogalvanic effect 
or the instantaneous value of the field F{t) in the CE- 
phase controlled AC current. This sensitivity to I{t) can 
be explored in a pump-probe setting, in which the in¬ 
stantaneous intensity can by manipulated by the delay 
between pump and probe pulses. The pump-probe de¬ 
lay may therefore serve as knob for fast charge trans¬ 
fer by the non-linear photogalvanic effect. To demon¬ 
strate this control we choose the intensity of both pump 
and probe pulse to be subcritical {Ii ^2 < ^c) with pump 
intensity /i = 2.4 x lO^^W/cm^ and probe intensity 
I 2 =0.6xl0^^W/cm^. However, the superimposed fields 
give rise to a maximum intensity of /max = 2.25/i = 
5.4x 10^^ W/cm^ above /c- The sign and amplitude of the 
induced current is controlled by the time delay between 
the laser pulses (Fig. 0]). For large positive and nega¬ 
tive delays, the transferred charge saturates at the same 
positive value. In contrast, for near-zero delay At = 0 
where the maximum intensity is attained, the DC current 
switches direction and the transferred charge becomes 
negative. Remarkably, during the period of strong over¬ 
lap the modulation of the DC current occurs for delays on 
the sub-fs time scale resulting from the strongly varying 
maximum instantaneous laser intensity as a function of 


pump-probe delay (Fig. IDd). Assuming that the charge 
transfer is governed by the central peak of the combined 
laser pulse, a simple estimate in analogy to Eq. [5] predicts 

D{M) = Sgn(Jc - Ji„ax(At)) ^nL , (7) 

where sgn denotes the sign function and /max (At) is 
the maximum instantaneous laser intensity for pump- 
probe delay At (Fig. IDd). In Eq. [3 we denote the non¬ 
linear generalization of the photogalvanic tensor by /3 nl • 
This simple model reproduces the temporal variation of 
D{At) in the full TD-DFT calculations remarkably well, 
underlining that the maximum instantaneous laser field 
drives the non-linear photogalvanic effect through tun¬ 
neling near the field maximum. 

In conclusion, we predict a non-linear extension of the 
photogalvanic effect into the strong-held regime giving 
rise to ultrafast DC currents in insulators illuminated by 
multi-femtosecond laser pulses. We observe a strongly 
non-linear intensity dependence and even a reversal of 
the induced currents above a critical intensity Ic asso¬ 
ciated with the transition from non-linear polarization 
currents to tunneling excitation. The charge transfer 
is rather insensitive to details of the laser pulse shape 
and carrier-envelope phase but strongly dependent on 
the maximum instantaneous held strength. The latter 
may be controlled by the pump-probe delay in a two- 
pulse setup giving rise to a distinct sign change in the 
transferred charge as function of the pump-probe delay. 
The non-linear photogalvanic ehect opens up opportuni¬ 
ties for light-held controlled femtosecond charge separa¬ 
tion with relatively modest requirements on the driving 
laser. Even many-cycle pulses without CEP stabiliza¬ 
tion can be used as the lattice structure instead of the 
CEP is employed to break the inversion symmetry along 
the laser polarization axis. The non-linear photogalvanic 
ehect is conceptually simpler than the CEP dependent 
charge transfer since no elaborate steering of the conduc¬ 
tion band electrons is necessary. Therefore, the ehect 
is robust against changes in the laser pulse parameters. 
This may be advantageous in particular for optical in¬ 
terconnects based on surface plasmon propagation [l8| 
where the pulse shape and duration of a surface plasmon 
wave packet is difficult to control. The sharp thresh¬ 
old intensity Ic for reversal of the current may provide 
a simple route towards femtosecond current switching 
and, moreover, a sensitive intensity calibration for laser 
pulses that directly measures the maximum electric held 
strength in the material. Finally, the photogalvanic ehect 
may also be investigated by associated terahertz emission 
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